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A luminous X-ray source is associated with a cluster (MGG-11) 
of young stars ~ 200 pc from the center of the starburst galaxy 
M82 1 2, The properties of the X-ray source are best explained 
by a black hole with a mass of at least 350 Mq 3, 4^ which is inter- 
mediate between stellar-mass and supermassive black holes. A 
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neeirby but somewhat more massive steir cluster (MGG-9) shows 
no evidence of such an intermediate mass black hole ^ 3^ raising 
the issue of just what physical characteristics of the clusters can 
account for this difference. Here we report numerical simula- 
tions of the evolution and the motions of steirs within the clus- 
ters, where steirs Eire allowed to mergers with each other. We 
find that for MGG-11 dynamical friction leads to the massive 
stars sinking rapidly to the center of the cluster to participate 
in a runaway collision, thereby producing a star of 800-3000 M©, 
which ultimately collapses to an black hole of intermediate mass. 
No such runaway occurs in the cluster MGG-9 because the leirger 
cluster radius leads to a mass-segregation timescale a factor of 
five longer than for MGG-11. 

Using the Keck NIRSPEC spectrometer, McCrady et al.^ have made accurate mea- 
surements of the bulk parameters of MGG-11 and MGG-9, two of the brightest star 
clusters in the central region of M82. Figure 1 shows X-ray^ and near-infrared^ images 
of the area of interest around the two clusters. The relative positional accuracies of 
both sets of observations are better than 1 arcsecond. However, the absolute pointing 
accuracy is much poorer, for both telescopes. Although apparently off-center, the po- 
sitions of the star cluster MGG-11 and the bright X-ray source are in fact consistent 
with one another (D. Pooley, private communication). 

Both MGG-9 and MGG-11 have quite similar ages, in the range 7-12 Myr i. The 
line-of-sight velocity dispersion of MGG-11 (cr^ ~ 11 A ± 0.8 km s~^) is somewhat 
smaller than that of MGG-9 {ur — 15.9 ± 0.8 km s~^). Combining these numbers 
with the projected half-light radii, r = 1.2 pc for MGG-11 and r = 2.6 pc for MGG-9, 
McCrady et al. estimate a total cluster mass of ~ (3.5 ± 0.7) x lO^M© for MGG-11, 
and about four times higher for MGG-9. 

In young dense clusters like MGG-11, supermassive stars may form through re- 
peated coUisions 6, 7, 8 xhe collision rate will be greatly enhanced if massive stars 
have time to reach the core before exploding as supernovae ^. Dynamical friction 
implies a characteristic time scale t^f for a massive star in a roughly circular orbit to 
sink from the half-mass radius R to the cluster center^O: 

(m) 0.138iV / R"^ V'"^ 

lOOM0ln(O.llM/lOOM0) \GM ) ' 

Here (m) and M are the mean stellar mass and the total mass of the cluster, respec- 
tively, is the number of stars, and G is the gravitational constant. For definiteness, 
we have evaluated tjf for a 100 Mq star. Less massive stars undergo weaker dynamical 
friction, and thus must start at smaller radii in order to reach the cluster center on a 
similar time scale. 

For MGG-11, we find tdf ~ 3 Myr, which is comparable to the main-sequence 
lifetimes of the most massive stars. On the other hand, for MGG-9, tdf ~ 15 Myr. 
Thus, massive stars in MCC-11 can easily reach the center of the cluster before 
exploding as supernovae, whereas those in MGG-9 do not. Given the high central 



2 



density of MGG-11, its massive stars, once accumulated in the cluster center, cannot 
avoid a runaway collision 6. 

We have tested this scenario by carrying out star-by-star simulations of MGG- 
11 and MGG-9 using two independently developed A^-body codes, Starlab^^' and 
NB0DY4 1'''' A more extensive discussion of those simulations is presented as Supple- 
mentary Information. Briefly, we choose the initial conditions of our model clusters so 
that today, at an age of 7-12 Myr, they have mass functions, luminosities, half-mass 
radii and velocity dispersions in agreement with the McCrady ct al. observations. 
Since we do not know either the initial or the current central densities of either clus- 
ter, the concentration c is treated as a free parameter controlling the initial central 
density of our model clusters. [The concentration parameter is defined as the ratio of 
the core radius to the tidal radius: c = log(i?t/i?c)-] 

We find that, for c > 2 (which for "King" is models is equivalent to a dimensionless 
central potential Wq ^ 9) our MGG-11 models do indeed show runaway growth via 
repeated collisions. Figure 2 presents, for several such simulations, the growth in 
mass of the star that would ultimately become the most massive star in the cluster. 
We will refer to this star simply as "the runaway star" . Based on detailed supernova 
calculations, we assume that stars having masses > 260 Mq collapse to black holes 
without significant mass loss in supernova explosions^^. Our stellar evolution models 
for stars with masses between 50 and 1000 M© are based on detailed calculations for 
such high-mass stars 

The quantitative differences evident in Figure 2 between the simulations performed 
with Starlab and those using NB0DY4 arc due mainly to the different radii assumed 
for the runaway star in those two packages. Stars of masses ^ 100 Mq in NB0DY4 are 
larger (in a time-averaged sense) by about a factor of 3 compared to those in Starlab. 
(More details are provided in Supplementary Information.) Because gravitational 
focusing dominates the collision cross section, this difference propagates linearly in 
the collision rate, explaining the factor ~ 3 difference in the final mass of the runaway 
star. Apart from this effect, we find that our qualitative results are quite insensitive 
to the details of the adopted evolution prescription. 

For MGG-11, the runaway star typically experienced a total of ~ 10-100 collisions. 
Most of them occurred during the first 3 Myr, that is, before the star became a black 
hole. The collision counterparts are usually 30-50 Mq main-sequence stars. By the 
time the runaway star collapsed to a black hole it had reached a mass of 800-3000 Mq. 
Later the black hole may capture a companion star to become an ultraluminous X-ray 
source ^9. 

No episode of runaway growth was seen in our MGG-11 models with c < 2, nor 
in any of the MGG-9 simulations, regardless of initial concentration. Thus we see 
very clearly that differences in bulk parameters, specifically the dynamical friction 
time scale for the most massive stars, can readily explain why MGG-11 might host 
an IMBH while MGG-9 does not. Simulations of MGG-11 with c > 2 (Wq ^ 9) 
reached core coUapse^o before IMyr. Models with c < 2 (Wq ^ 8) also showed some 
increase in central density, but the maximum density was much lower (and occurred 
later) because the collapse was stopped by supernova mass loss after ~ 3 Myr. Thus, 
a collision runaway could not occur in these clusters. 
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Figures summarizes the results of our simulations, illustrating how both high 
initial concentration and short dynamical friction time scales are needed to lead to a 
collision runaway. In addition to the parameters for MGG-9 and MGG-11, we also 
indicate on the figure the location of the young star cluster [W99]Pi, one of several 
star clusters in the Antennae system for which accurate structure parameters have 
been determined ^. Its parameters are quite similar to those of MGG-9. The cluster 
[W99]l has a bright i^o.2-iokcV — lO"'^'^ erg/s X-ray point source as a counterpart, ^3 the 
luminosity of which is consistent with a "normal" high-mass X-ray binary containing a 
magnetized neutron star or a stellar-mass black hole ^3. Star clusters such as MGG-11, 
MGG-9 and [W99]l are richly populated with black holes and neutron stars — based 
on our adopted mass function, we expect these clusters to contain some 1400 stellar- 
mass black holes and up to about 1000 neutron stars. The formation of an ordinary 
high-mass X-ray binary is therefore not surprising. 

The relevant parameters of five other young clusters in the Antennae^ all fall far 
to the right of Figure 3 and are therefore not expected to contain IMBHs. This is 
consistent with the absence of ultr aluminous X-ray sources in these clusters. The 
Milky Way galaxy contains (at least) four young dense star clusters, of which the 
Arches, Quintuplet, and NGC 3603 have the right conditions for multiple stellar 
collisions to occur. However, the relatively small number of stars in these systems 
may prevent the growth of an object massive enough to collapse to an IMBH. In 
any case, these clusters are currently too young to have experienced any supernovae. 
The slightly older star cluster Westerlund 1^^ is sufficiently massive, and fulfills the 
criteria for producing an IMBH via the process described above. So far, no bright 
X-ray source has been found in this cluster. 

Based on theoretical considerations, supported by extensive numerical simulations, 
we have shown that an IMBH of mass ~800-3000 Mq is expected to form through 
runaway collisions in the star cluster MGG-11, but not in MGG-9. The requirement 
for the formation of such an IMBH in MGG-11 is that the cluster was born with c ^ 2 
and tdi ^ 4Myr. While high by the standards of typical open and globular clusters 
observed in our Galaxy, such a density is not uncommon among young ( ^ 3 Myr old) 
star clusters. Examples include NGC 3603 (c > 2.08)25 in the Milky Way Galaxy, 
and R136 in the 30 Doradus region of the LMC ^6. We propose that the IMBH thus 
formed is the origin of the LQ 2-iokeV — 10^^ erg/s ultraluminous X-ray point source 
seen at the position of MGG-11. 
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Figure 1: ChandraX-raj image (in R.A. and Dec.) of the relevant part of the starburst 
galaxy M82 with the observed star clusters indicated. The color image is from the 28 
October 1999 X-ray observation by Matsumoto et al. (2001), 3 with about arcsecond 
positional accuracy. The brightest X-ray source (M82 X-1) is near the center of the 
image. The star clusters, from table 3 of McCrady et al. (2003), ^ are indicated by 
circles. The positions of the two star clusters MGG-9 and MGG-11 are indicated with 
squares. The magnified infrared images of these star clusters from the McCrady et 
al observations are presented in the upper right (MGG-11) and lower left (MGG-9) 
corners. A recently discovered 54.4 ± 0.9 mHz quasi-periodic oscillator is not shown 
because of its low (7 arcsecond) positional accuracy^, but its position is consistent 
with the X-ray source in MGG-11. A millimeter source^''' roughly centered around 
the two clusters is not shown either, because it is a large shell-like structure with a 
diameter of ~ 14" x 9". 
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t [Myr] 

Figure 2: The growth in mass of the colhsion runaway star with time. The choice 
of initial concentration is labeled by the central potential, where W12(9) implies 
Wq — 12(9). For both choices, the top curves give the NB0DY4 results, and the 
bottom curves the Starlab results. The runaway masses in NB0DY4 are larger, since 
that code adopts larger stellar radii, as discussed in the Supplementary material. 
The star symbols indicate the moment when the runaway experiences a supernova, 
typically around 3 Myr. The open and filled stars indicate simulations performed 
with NB0DY4 and Starlab, respectively. The sohd and dashed curves show for a 
Salpeter2 IMF with a lower limit of 1 Mq and c ~ 2.1 {Wo = 9) and c ~ 2.7 {Wq = 12). 
The dash-dotted curves are for two models with Wq = 9 with an upper limit to the 
IMF of 50 M0, instead of the standard 100 M0 used in the other calculations; we 
terminated these runs at the moment the runaway star experiences a supernova. The 
dash-3-dotted curve shows the result for Wq — 12 with a Salpeter IMF and with 
10% primordial binaries. Finally, the dotted curve shows results for Wq — 9 and a 
Kroupa26 IMF with a minimum mass of 0.1 M©, in a simulation with 585,000 stars. 
The observed age range of MGG-11 and MGG-9 is indicated by the horizontal bar 
near the bottom of the figure. 
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Figure 3: The area of parameter space for which runaway colhsion can occur, and 
where the process is prevented. Conditions for runaway merging identified through 
our simulations in the {tdf, c} plane, where t^f is the dynamical friction time scale for 
100 M0 stars, and c the cluster concentration parameter. The horizontal error bars 
for the initial conditions for MGG-9 and MGG-11 reflect errors in the observed cluster 
mass and projected half light radius ^. The best fit for the concentration parameter of 
the Antennae star cluster [W99]l is indicated by a vertical bar near t^f — 17Myr, the 
horizontal error bar reflect the uncertainty in the measured cluster mass and radius 
9. Solid circles indicate simulations resulting in runaway merging, while open circles 
correspond to simulations in which no runaway merging occurred. The evolution of 
the mass of the runaway for the two leftmost flUed circles are presented in Figure 
2, for a variety of initial conditions. Star clusters in the upper left corner enclosed 
by the dashed line are expected to host an IMBH formed by the runaway growth of 
a single star. In the other parts of the diagram (lower concentration and/or larger 
dynamical friction time) an IMBH cannot form by this process. 
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Supplementary material 



This part of the paper will be published online and it has a 
separate bibliography 

In this section we discuss the selection of initial conditions, numerical methods, 
and the results of our simulations of the star clusters MGG-9 and MGG-11. 

The first step in modeling this problem is to reconstruct the conditions under 
which MGG-11 was born 7-12 Myr ago. McGrady et aP determine the cluster age 
and mass-to-light ratio, and from these conclude that the mass function is consistent 
with a power- law having the Salpeter^ slope and a lower mass limit of ~ 1 Mq ^. At the 
current age of the cluster, all stars more massive than ~17-25Mq3 have experienced 
supernovae. We adopt this mass function in most of our simulations, with an upper 
limit of 100 M0. This choice of mass limit is motivated by the fact that there are 
several clusters in the Milky Way Galaxy with characteristics similar to the two 
clusters studied here (MGG-11 and MGG-9) that are young enough for the highest- 
mass stars still to be present. These clusters are the Arches^, Quintuplet^, NGC36036 
and the central cluster R136'^ in the 30 Doradus region in the large Magellanic cloud 
(LMC). The publications cited here support the existence of stars of ^ 100 Mq in 
these clusters. In addition, an upper mass limit of ~ 100 Mq in the Galactic disc is 
probably quite realistic ^. 

The 1 Mq lower mass cut-off for the initial mass function (IMF) is somewhat 
controversial. However, the mass functions in several of the young star clusters in the 
Antennae galaxies^, the Arches cluster^. Quintuplet clustcrio, NGC1705-lii and the 
star cluster MGG-F in M82ii are consistent with a similar lower mass cut-off. For 
our study we adopt the observed mass function as a first choice, but in later runs we 
vary the mass function to test the effect of changing it. 

At zero age the mean mass of our adopted IMF is (m) ~ 3.1 Mq. At the observed 
age of 7-12 Myr the mean stellar mass has dropped to (m) ~ 2.5 — 2.7 Mq 3. What 
actually happens to stars of ~17-25 Mq which experience a supernova before 7-12 Myr 
is not crucial for this estimate, as there are only 1500-2500 such stars. If we naively 
assume that all the mass in these stars is simply lost from the cluster we obtain mean 
stellar mass of 2.5-2.6 Mq. If we assume that stars more massive than 25 Mq turn 
into 10 Mq black holes and stars of lower mass into Chandrasckhar [1.4 Mq] mass 
neutron stars, the mean stellar mass drops to (m) ~ 2.6 — 2.7Mq. With today's 
measured mass of M = 3.5 x 10^ MqI the cluster then contains N = 130, 000-140,000 
stars, implying an initial total mass of Mq — 4.1-4.5 x 10^ Mq. 

In order to model the cluster we must know its initial half-mass radius R. We 
obtain this parameter by de-projecting the observed half-light radius and correcting 
for the expansion of the cluster during its lifetime. The two corrections are compa- 
rable, but with opposite sign: the half-mass radius is ~25 percent larger than the 
observed half-hght radius, and mass loss due to stellar evolution drives adiabatic ex- 
pansion of the cluster radius by ~20 percent. The initial half-mass radius of MGG-11 
is then R ~ 1.2 pc. Here we assume that mass follows light after 12 Myr, i.e. that the 
light profile is not significantly affected by the segregation of the most massive stars. 
This turns out to be a reasonable assumption. Note that the mass-to-hght ratio for 
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MGG-11 reported by McCrady et al. is L/M = 3.5t?:5, and that adding a 3000 
non-luminous object to the cluster increases this value only by about 1 percent. A 
similar analysis for MGG-9 results in a total mass at birth of Mq ~ 1.8 x 10^ Mq and 
an initial half-mass radius of i? ~ 2.6 pc. 

We ignore the effect of the tidal field of M82 in our simulations. The runaway 
growth studied for these clusters takes place within the dynamical friction time scale 
of the most massive stars and the process discussed here takes place deep within 
the cluster core. Even though the cluster may lose a few stars to its host Galaxy, 
this does not affect the processes driving the collision runaway. Note that the two 
clusters MGG-9 and MGG-11 are close together on the sky, and could in principle be 
bound in a binary cluster. At the 3.9 Mpc^^ distance of M82, the projected distance 
between MGG-9 and MGG-11 is about 30 pc. We note that in the LMC there arc a 
few hundred binary clusters for which the projected separation is smaller than 20 pc 

12, 13 

Our simulations were performed using two independently developed A^-body pack- 
ages, Starlab^^' q^^^ NB0DY4 Both codes achieve their highest speed when 
running on the GRAPE-6 special-purpose hardware ^o. A comparison between 
Starlab and NB0DY4 is presented by Hcggic (2001)2i. 

Both codes incorporate parametrized stellar evolution and allow for the possibility 
of direct physical collisions between stars, thus including the two key physical elements 
in our runaway scenario. The criteria for stellar collisions in both codes are quite 
similar; a collision occurs if the distance between two stars becomes smaller than 
the sum of the stellar radii. During a collision mass and momentum are conserved. 
These arc reasonable assumptions since the relative velocity of any two colliding stars 
is typically much smaller than the escape speed from either stellar surface 22. As the 
runaway star grows in mass its angular momentum increases due to off-axis collisions. 
Our stellar evolution prescriptions did not incorporate special treatments to account 
for this effect. Rapidly rotating stars tend to develop significant horizontal turbulence 
which may dramatically affect the further evolution of the star 23, 24 Galculations for 
a rapidly rotating star of 120 M0 results in a ~ 15 percent extension of the lifetime 
compared with a non- rotating star 25. This effect is stronger for lower mass stars, but 
has not yet been explored for more massive stars. 

The collision radius of a star is taken to be identical to the stellar size, except 
for black holes for which we use the tidal radius, which is considerably larger than 
the Schwarzschild radius. The zero-age radii of various stellar evolution models are 
presented in Figure 5. The radii in Starlab are systematically smaller than those in 
NB0DY4. The collision cross sections in the appropriate range of initial conditions are 
dominated by gravitational focusing, and therefore the stellar size propagates linearly 
in the collision rate, and consequently also in the final mass of the runaway star. 

We simulated MGG-11 and MGG-9 using initial mass functions and half mass 
radii based on the observations of McGrady et al.^ Most of these calculations were 
performed with 131072 (128k) stars. The initial density distribution is poorly con- 
strained by the observations and we therefore performed a series of simulations with 
various choices for the initial density profile, ranging from the very shallow c ~ 0.67 
(equivalent to Wq = 3) to highly concentrated c ~ 3.4 (equivalent to Wq = 15). To 
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further test our theoretical understanding that a coUision runaway is prevented by 
the first supernova explosions for t^f ~ 4 Myr we systematically increased the initial 
tdf for two selected values of c in our MGG-11 simulations (by varying the initial 
half-mass radius R), until the collision runaway was prevented. The results are pre- 
sented in Figure 3 of the main text. Note that we plot the horizontal error bars for 
both clusters near the top of the figure for practical reasons; these should not be 
interpreted as upper limits on the concentration parameter. 

In order to illustrate the early core collapse in the more concentrated models, 
and the absence thereof in the shallower models, Figure 4 plots the evolution of 
the core radius for a selected number of simulations for MGG-11. We discuss each 
curve in turn, starting at the top. The core radius of the shallow model, c ~ 0.67 
(Wo = 3), hardly changes with time. The intermediate model c ~ 1.8 (Wq = 8) 
almost experiences core collapse near t = 3 Myr but as stellar mass loss starts to 
drive the expansion of the core it never really experiences collapse. Inspection of 
Figure 3 of the main text shows that neither of these clusters experienced a collision 
runaway. Gore collapse occurs in the c ~ 2.1 (Wq — 9) model near t — 0.8 Myr, 
which also corresponds to the time of the first collision. The c ~ 2.7 (Wq = 12) 
simulation is so concentrated that it starts virtually in core collapse, and the entire 
cluster evolution is dominated by a post-collapse phase. Inspection of Figure 2 in the 
main text reveals that in these calculations the runaway also tends to start earlier. 
The larger fluctuations in the core radius for models Wq — 9 and Wq — 12 compared 
to the less concentrated models reflect the orbit of the runaway star around the 
density center of the star cluster. The runaway star is always found in the core, but 
occasionally at a substantial distance from the center. 

We performed four simulations for MGG-11 adopting Wq ~ 9 with different initial 
mass functions in order to study the effect of the mass function on the calculations. 
Reducing the upper limit on the IMF from 100 M0 to 50 M0 reduces the mass of 
the runaway star at the moment of the supernova explosion by about 20 percent. 
Simulations with the Kroupa (2001)26 IMF and a lower mass limit of 0.1 M©, as well 
as one run with a Salpeter IMF and a lower mass hmit of 0.2 Mq, resulted in a ~ 30% 
lower runaway mass. These latter simulations were run with 585,000 stars. The lower 
mass of the runaway star in the simulations with the reduced upper mass limit is a 
direct consequence of the smaller number of massive stars in the simulation, causing 
the average mass of the collision counterparts to be smaller and the collision rate to 
be lower. We also performed a calculation using the initial conditions for MGG-11 
but with the zero-age stellar radii pubhshed by Ishii et al. (1999)^3 for the collision 
runaway. This calculation resulted in the a similar black hole mass as in the NB0DY4 
calculations with otherwise the same initial conditions. 

In addition, we performed several simulations for MGG-9, which is well outside 
the regime where a collision runaway can occur (see Figure 3 of the main text). 
For MGG-9 the IMF is consistent with a Salpeter IMF with a lower mass limit of 
0.2 M0, and also with the Kroupa IMF. However, due to computational limitations we 
could only perform simulations with the 1 Mq cut-off IMF for MGG-9. Simulations 
with the Kroupa IMF would require 3.6 million stars, which is beyond our current 
capabiUties. In the case of the Salpeter IMF with a lower limit of 1 Mq, and with other 
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parameters chosen consistent with the observations of McCrady et al for MGG-9, we 
did not observe runaway merging. 

Finally, we performed two calculations for MGG-11 using c ~ 2.7 {Wq = 12), and 
including 10 percent primordial binaries. The current binary fractions for MGG-11 is 
not determined by the observations, but based on studies of other clusters, 10 percent 
seems a reasonable choice 27, 28, 29 select initial orbital separations ranging from 
Roche- lobe contact (at a binding energy of about 300 kT) to that corresponding to 
a binding energy of ~ 3kT, which is near the dividing line between hard and soft 
binaries ^o. (Note that the total kinetic energy of the stellar system is ^NkT.) The 
other initial binary parameters were selected as in Portegies Zwart et al. (2001) 
15. Though modest, this fraction should be sufficient to completely dominate the 
dynamical evolution of the cluster core ^i. 

The presence of primordial binaries in our models causes an increase in the overall 
merger rate by about a factor of four. Most of these mergers do not involve the run- 
away star but occur as a result of unstable mass transfer in a primordial binary. The 
collision rate with the runaway star increases by about 30 percent, and as a conse- 
quence the mass of the runaway star at the time of supernova is 22-36 percent larger 
than without primordial binaries. We have neglected the possibility of primordial 
mass segregation, which would further favor the runaway process, as massive stars 
would be born preferentially in the cluster core. 
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Figure 4: Evolution of the core radius for four simulations of the star cluster MGG- 
11. These calculations are performed with Starlab with c ~ 0.67 {Wq = 3), c ~ 1.8 
(H/q = 8), c ~ 2.1 {Wo = 9, bold curve) and c ~ 2.7 {Wq = 12, dotted curve) 
indicated along the right edge of the figure as W3, WS, W9, and W12 respectively. 
The W9 curve is plotted with a heavy line to distinguish it from the curves for W8 
and W12. The age range of the observed clusters is indicated near the bottom of the 
figure. The large variation in the core radius for the Wq = 9 and Wq = 12 models is 
caused in part by the presence of the runaway star. 
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Figure 5: Stellar radius as a function of zero-age mass for various models of high- 
mass stars. The lower and upper solid lines indicate the radii of zero-age and terminal 

age main-sequence stars from the stellar evolution model used in Starlab, which are 
based on the prescription of Eggleton ct al. (1989)3. The dashed lines are taken from 
the stellar evolution prescription used in NB0DY4 which are based on the evolutionary 
calculation by Hurley et al. (2000)32. The triangles are from the calculations of Ishii 
et al. (1999)33, the squares are taken from Stothers (1997) 34. 
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